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Abstract. We describe the results and the current status of our program to 
search for intracluster planetary nebulae (ICPNe) in the Virgo and the Coma 
clusters of galaxies. For the Virgo cluster, ICPNe candidates were detected with 
narrow band imaging combined with broad band imaging. The use of two narrow 
band filters for Ha and [O III] lines tuned for the cluster redshift enabled us to 
detect secure ICPNe candidates. Kinematics of the detected ICPNe confirm that 
the Virgo cluster is in a highly unrelaxed state. In order to detect ICPNe in the 
Coma cluster, which is more than five times more distant than the Virgo cluster, 
we devised a multi-slit imaging spectroscopy technique and successfully detected 
35 secure ICPNe candidates for the first time in this cluster. We have found a 
hint that the Coma cluster is currently in the midst of a subcluster merger. 



Introduction 



Cosmological simulations predict that a substantial fraction of stars in galaxies 
were stripped off during the assembly of a cluster of galaxies. These stars were 
dispersed into intracluster space and now recognized as the intracluster stellar 
population. This population is a sensitive probe of cluster assembly, and more 
generally, structure formation of the universe. The total amount, spatial distri- 
bution (morphology), kinematics, and metallicty are among the key parameters 
which characterize the population. 

The intracluster ste l lar population has been observed as diffuse intracluster 
light dMihos et af]l2005l: iKrick et al.ll2006lh or as individ ual stars, i.e., plane- 
tary nebulae (Arnaboldi et al. 1996 ; lOkamura et al.l 2002), or as red giant stars 
( Ferguson et al. 19981 ). The intracluster planetary nebulae (ICPNe) are the only 
probe that enables us to investigate the kinematics of the intracluster stellar 
population. We launched a program to search for ICPNe in clusters of galaxies 
with a range of properties on the advent of Subaru telescope and its wide-field 
Camera, Suprime-Cam. 



2. Observing Strategy 

Initial targets of our program are the Virgo cluster at the distance of (to — M)q ~ 
31 and the Coma cluster at (m — M)q ~ 35. The Virgo cluster is irregular and 
loose in shape, and dynamically unrelaxed containing a considerable fraction 
of spiral galaxies, while the Coma cluster is regular and rich in shape, and 
dynamically more relaxed than the Virgo cluster. 
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We applied a conventional narrow band plus broad band imaging technique 
to the nearby Virgo cluster. Unique feature of our search is the use of two narrow 
band filters for Ha and [O III] lines tuned for the cluster redsh ift. Previous 



plane tary nebula searches over the image of Virgo ellipticals (e.g, iJacobv et al 



1990, hereafter JCF) used the strongest [O III] line only. The PNe candidates 
selected by the [O III] single line technique is contaminated by foreground and 
background emission line objects. Especially, when the single line technique is 
used to search for ICPNe in the 'blank region' free from Virgo member galaxies, 
the contamination could dominate the candidates. Use of the two lines and 
broad band images enabled us to identify secure ICPNe candidates. However, 
the follow-up spectroscopy is necessary to the final confirmation and kinematic 
studies. We used Telescopio Nazionale Galileo (TNG) 4 m telescope and Very 
Large Telescope (VLT) for the follow-up spectroscopy. 

No imaging technique can be successfully applied to the detection of ICPNe 
in the distant Coma cluster. At the Coma distance, even the brightest PN is 
buried in the sky noise and undetectable in the narrow band image. Increase of 
the signal-to-noise ratio of the emission lines by reducing the sky noise is crit- 
ically important. To achieve this, we devised a multi-slit imaging spectroscopy 
(MSIS) technique, which was sucessfully applied to the Coma cluster. 



3. Virgo Cluster 

3.1. Imaging Observation and Selection Criteria 

Observations, sel ection criteria for secu re P Ne candidates, and in itial results 
were described in lOkamura et all ( 20021 ) and lArnaboldi et al.l ( 20031 ). Only the 



outline is repeated here. 

In 2001 March - April a field in the Virgo Cluster at (a, 5)(J2000) = 
(12 /l 25 m 47. s 0, +12°43 / 58"), just sou th of M84-M86 was observed during the 
commissioning of the Suprime- Cam (Miyazaki et al.ll2002l ). at the prime focus 



of the Subaru 8.2 m telescope ( Ive et al. 20041 ). Suprime-Cam covers an area 



of 34' x 27' of the sky, with a resolution of 0."2 pixel 1 . The field was imaged 
through two narrow band and two broadband (V and R) filters. The two narrow 
band filters have (A c , AA) = (502lA, 74A) and (6607A, lOlA), corresponding to 
the [O III] and Ha emissions at the redshift of the Virgo cluster. Exposure times 
are 900, 720, 3600, and 8728 seconds for V,R, [O III], and Ha, respectively. 

We use the two-color diagram [O III] - Ha versus [O III] - (V+R) to select 
our emission-line objects. Simulated populations of PNe, single-lined emitters, 
and continuum sources were constructed and used to outline the regions in the 
two-color diagram inhabited by the different kinds of objects (Fig. 1, left panel). 
On the basis of this simulation, we set selection criteria. Two bright galaxies 
M86 (NGC 4406) and M84 (NGC 4374) appear in our field. These two galaxies 
were surveyed for PNe by JCF, and we can evaluate our selection criteria against 
the JCF photometry of matched samples. The right panel of Fig. 1 shows the 
comparison for M84 field, where 74 objects are matched. We found that at least 
% contamination is present in the JCF sample selected by [OIII] line only. 



Intracluster Planetary Nebulae 



3 




secure PN candidates 



. / 



dpntinuum objects 

/ 



lik ely PN candidates 



Compact HII regions 



[OIII] - (V+R) 



plus galaxies 



JOIII] - (V+B) 



Figure 1. (left) Two-color diagram of the simulated population of point 
sources. Dots are PNe, stars are single emission-line objects, and triangles 
are continuum objects, (right) Two-color diagram for the M84 matched sam- 
ple. Dots are [O III] sources matched to the JCF sample. Triangles (19 
objects) indicate the best PN candidates, and circles (26 objects) indicate the 
single-line emitters, most likely high excitation PNe in this situation. Squares 
indicate continuum objects. Three pentagons indicate objects with unreliable 
colors, and one of them was fou nd to be a compact H II region apparently 
floating in the intracluster spa ce dGerhard et al.|[2 002'l . Our selection criteria 
are shown by the thick lines. ([Arnaboldi et al.ll2003f ) 



3.2. Follow-up Spectroscopy for the Virgo Cluster 



We carried out follow-up spectroscopy of t he ICPN candidates w ith multifiber 
FLAMES spectrograph at UT2 on VLT (|Arnaboldi et alJl2004l ). Candidates 



were drawn from the three survey fields, FCJ, CORE, and SUB, as shown in 
the left panel of Fig. 2. Spectra covered a wavelength range of 500 A, centered 
on 4797A with a resolution of A/AA = 7500. This gives a velocity resolution of 
40 km s~ x and a typical velocity error of 12 km s^ 1 . 

A large fraction of the photometric candidates with m(5007)<27.2 were 
confirmed and A4959 line of the [O III] doublet was detected for the first time 
in about a half of confirmed ICPNe. A total of 40 ICPNe candidates (15/12/13 
for FCJ/CORE/SUB fields) were confirmed, and their radial velocities were 
measured. 

Radial velocity distributions of the ICPNe in the three fields are shown 
in the right panel of Fig. 2. They are considerably different from each other. 
This confirms the view that Virgo is a highly nonuniform and unrelaxed galaxy 
cluster, consisting of several subunits that have not yet had time to come to 
equilibrium in a common gravitational potential. 
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Figure 2. (left) Three fields in the Virgo cluster where follow-up spec- 
troscopy was made. Circles show FLAMES field of view, (right) Rad ial 
velocity distributions of ICPNe in the three fields. (jArnaboldi et al.ll2004f ) 



Toward the Coma Cluster 



4.1. MSIS Technique 



As already mentioned, no imaging technique can be successfully applied to the 
detection of PNe in the Coma cluster. We devised a technique that is simi- 
lar to the approach used to search for Lya-emitting galaxies at very high red- 
shift ([Stern Sz Spinrad 1999). It combines multi-slit grism spectroscopy with 



a narrow band filter centered around the [O III] A5007 line at the redshift of 
the C oma cluster. We use the Faint Object Camera and Spectrograph (FO- 
CAS; iKashikawa et al.l 120021 ) at the Cassegrain focus of Subaru telescope (see 



Gerhard et al.l 120051 for the method and initial results). 

A mask of parallel multiple slits is used to obtain spectra of all PNe that 
happen to lie behind the slits. Because the [O III] emission lines from PNs 
are only a few kilometers per second wide, their entire flux still falls into a 
small number of pixels in the two-dimensional spectrum, determined by the slit 
width and seeing. On the contrary, the sky emission is dispersed in wavelength, 
allowing a large increase in the signal-to-noise ratio (S/N). The narrow band 
filter limits the length of the spectra on the CCD detector, so that many slits 
can be exposed. This technique, which is a blind search, is referred to as the 
multi-slit imaging spectroscopy (MSIS) technique. No conventional imaging 
technique can decrease the sky surface brightness in a similar way. 

Figure 3 shows the MSIS field and the slit mask. The circular FOCAS field 
of view (FOV) has 6' in diameter (166 kpc at the assumed distance of 95 Mpc). 
The spatial resolution of FOCAS is 0."1 pixel -1 , so the 6' diameter corresponds 
to 3600 pixels. The light passing through the narrow band filter and a 0."6 wide 
long slit, and dispersed by the grism, projects down to a spectrum of about 43 
pixels. A mask was therefore manufactured with uniform long slits spaced every 
50 pixels and interrupted only by short sections to ensure mechanical stability. 
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Figure 3. MSIS techniqu e, (left) Location of the MSIS field, (right) Slit 
mask. (jGerhard et alJl2005f ) 



We used 300B grating, which gives a measured dispersion of 1.45A pixel 1 on 
the two FOCAS CCD chips. The effective spectral resolution is 7.3A, or 440 
km s -1 . We also used the N512 filter with FWHM 60A, centered at A c = 
512lA, the wavelength of the redshifted [O III] emis sion from a PN at the mean 
recession velocity of the Coma Cluster (6853 km s~ 1 : IColless h Dunnlll996l ). The 
60A FWHM includes only ±1.6cr the velocity dispersion of galaxies. By taking 
exposures at 8 positions sucessively shifted along the dispersion direction, we 
can roughly cover the whole field of view. 



4.2. Hint of On-going Subcluster Merger in the Coma Cluster 

Observations at three mask positions centered on the field at (a, 6) (J2000) = 
(12 fe 59 m 41. s 8, +27 53'25."4) (see Fig. 3) were carried out during the nights of 
2004 April 21-23 . We detected a total of 60 emission line objects. Among them, 
35 are ICPNe candidates with no continuum flux, 5 are PN candidates projected 
against the extended halos of Coma member galaxies, and 20 are background 
galaxies. Details of the data re duction and a catalog of ICPNe candidates are 
given m lArnaboldi etaD (|2007l ). 



Figure 4 shows the velocity histograms of ICPNe (left), Coma galaxies in the 
MSIS field (middle), and background emission- line galaxies (right). We detect 
clear velocity substructures of ICPNe, i.e., intracluster stellar population, within 
a 6 arcmin diameter MSIS field. A substructure is present at ~ 5000 km s _1 , 
which is probably from infall of a galaxy group, while the main component is 
centered around ~ 6500 km s _1 , ~ 700 km s _1 offset from the nearby cD galaxy 
NGC 4874. The kinemati cs of ICPNe found in this s tudy and morphology of the 
diffuse intracluster light ( Thuan Kormendv 19971 ) show that the cluster core 



is in a highly dynamically evolving state. In combination with galaxy redshift 
and X-ray data this argues strongly that the cluster is currently in the midst of 
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Figure 4. Velocity histograms of ICPNe (left), Coma galaxies in the MSIS 
field (middle) , and background emission- line galaxies (right) . Filter transmis- 
sion curve is shown in the left and right panels. (Gerhard et al. 2007) 



a subcluster merger during which the elongated distribution of diffuse light has 
been created. The two subcluster cores (NGC 4874 and NGC 4889) are presently 
just after their second cl ose passage. The lik ely orbits of the subcluster cores 
are illustrated in Fig. 5. ( Gerhard et al.l [20071 ) 




Figure 5. The likely orbits of NGC 4889 and NGC 4874 up to their present 
positions are sketched. Contours sh ow the surface bright n ess d istribution of 
diffuse intracluster light t aken from iThuan fc Kormendvl (|1997l ). The circle 
indicates the MSIS field. (|Gerhard et al.ll2007t) 
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4.3. New Deep Imaging Program for Diffuse Intracluster Light 

Inspired by the work bv lMihos et~aD (120051 ) and lKrick et~aD (|2006l ). we started 
a new very deep imaging program to detect the diffuse intracluster light in a 
wide 50' x 50' field centered on the Coma cluster using a CCD camera attached 
to the 105-cm Schmidt telescope at Kiso observatory, Institute of Astronomy, 
University of Tokyo. We hope to obtain some data in spring, 2008. 
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Questions and Answers 



M. Strauss: How large would the contamination of the Virgo ICPNe sample 
be by background galaxies if you didn't have the Ha filter? 

S. Okamura: Lower limits of the fraction of contamination estimated from 
the comparison with FCJ are 39% and 56% for M84 sample and M86 sample, 
respectively. In a blank intracluster field, the contamination could be much 
larger. So, the use of Ha filter is critically important. 

N. Kawai: Can you tell whether these ICPNe were originally formed in the 
member galaxies or formed outside? 

S. Okamura: Most of them are likely to be formed in the member galaxies 
and stripped off later during the cluster assembly. It is, however, still an open 
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question whether or not some ICPNe were formed outside galaxies in the very 
early phase of cluster formation. 

J. Krick: G. Bernstein et al. measured the amount of diffuse intracluster 
light in the Coma cluster to be ~ 50% of the total cluster light. How do your 
values from the ICPNe compare to this? 

S. Okamura: We have not yet made the estimate for the Coma cluster. It 
is probably highly unreliable if we make the estimate based on the number of 
ICPNe detected in the present study over a fraction (~3/8) of a single small 
MSIS field. I hope our new deep imaging program with Kiso Schmidt telescope 
would give a reliable estimate. 

N. Tamura: Any estimation of metallicity for ICPNe? 

S. Okamura: Estimation of metallicity is a very difficult job yet to be done. 



